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NOTATION 


Wave  frequency,  cycles  *  sec’^  or  Hertz 

Frequency  corresponding  to  the  peak  of  the  wave  spectrum, 
cycles  •  sec“l 

•2 

Acceleration  due  to  gravity,  9.8087  m  •  sec 
Spectral  moment  of  order  zero 
Long-crested  wave  spectral  density  ordinates 
Wind  speed,  at  10  m  above  ocean  surface 
Fetch 

Significant  wave  height,  average  of  one-third  highest 
double  amplitudes 

Modal  wave  period,  period  corresponding  to  the  frequency 
of  the  peak  of  wave  spectrum 

Phi  1 1  ip's  constant 

A  constant  dependent  on  the  significant  wave  height  and 
modal  wave  period 

Ratio  of  the  maximal  spectral  energy  to  the  maximum  of  the 
corresponding  Pierson-Moskowi tz  spectrum 

Left,  right  side  widths  of  spectral  peak  of  JONSWAP 
spectra 
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ABSTRACT 


For  simplicity  as  well  as  consistency  with  the 
current  state-of-the-art  in  seakeeping  performance 
assessment,  a  wave  spectral  formulation  for  fetch- 
limited  ocean  areas  which  is  dependent  only  on  the 
significant  wave  height  and  modal  wave  period  Is 
desirable.  A  modified  version  of  the  JONSWAP 
spectrum  is  therefore  derived,  but  as  is  the  case 
with  the  usual  fetch-dependent  JONSWAP  spectrum,  it 
contains  too  much  energy  for  fetches  above  about  40 
nautical  miles.  This  inconsistency  is  probably  due 
to  certain  parameter  relationships  which  are  all 
based  on  least  squares  fits.  Therefore,  a  new 
parameter,  6,  is  developed  to  replace  the  usual  a 
parameter,  correcting  for  the  parameter's  nonuni- 
versal ity. 
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Number  SF  43  421  001.  It  is  identified  by  Work  Unit  Number  1504-100  and 
1500-104,  respectively,  at  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC). 


INTRODUCTION 


The  primary  goals  of  the  Joint  North  Sea  Wave  Project  (JONSWAP),  which 
originated  in  1967,  were  to  measure  the  growth  of  waves  under  limited  fetch 
conditions  and  to  analyze  attenuation  of  waves  propagating  into  shallow 
water,  see  References  1  and  2.*  The  fetch  dependence  of  the  measured  one¬ 
dimensional  frequency  spectra  was  investigated  by  parameterizing  them  with 
an  analytic  function  derived  by  least  squares  fit  techniques.  The  resulting 
function,  currently  the  most  widely  used  spectrum  for  representing  fetch- 
limited  seas,  is  known  as  the  JONSWAP  spectral  density  equation  and  is 
given  by  _ 
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which  reflects  the  five  parameters  of  fQ  a,  y,  and  shown  in  Figure 

1.  f  is  the  wave  frequency  in  cycles  •  sec  * ,  f  Is  the  frequency  at  the 
spectral  peak,  and  g  is  the  acceleration  due  to  gravity.  In  this  report, 
a  "mean"  JONSWAP  spectrum  has  been  used,  so  that  y  is  3 - 3 »  <?_  is  0.07  and 
Ojj  Is  0.09.  The  two  parameters,  a  and  f  ,  are  dependent  on  the  wind  speed 
and  fetch  such  that 

a  =  0.076  x‘°*22  (2) 


and 


(3) 


where 


and 


fQ  -  3.5  X"0,33,  X  <  104 


W 

(5) 


X  is  the  fetch  in  nautical  miles.  The  wind  speed  U  is  taken  to  be  at  an 
elevation  of  10  m  above  the  surface  and  Is  in  units  of  knots. 

The  JONSWAP  spectral  form  represents  a  generalization  of  the  Pierson- 
Moskowitz  form  by  inclusion  of  fetch  as  an  additional  parameter  to  wind 
speed.  If  o  is  0.0081  and  y  is  1  in  equation  (1),  the  JONSWAP  spectral 
form  is  identical  to  the  Pierson-Moskowi tz  form.  In  general,  the  JONSWAP 
spectrum  contains  more  peak  energy  than  the  corresponding  Pierson-Moskowi tz 
spectrum  for  the  same  values  of  a  and  f  ,  see  Figure  1. 


MODIFICATION  OF  THE  JONSWAP  SPECTRUM 
In  previously  completed  ship  seakeeping  analyses,  e.g.,  see  Reference 
3,  in  which  the  ship's  operations  in  fetch-limited  waters  has  to  be 
assessed,  it  has  been  customary  to  apply  the  JONSWAP  spectrum,  defined  by 
wind  speed  and  fetch  as  given  in  equation  (1).  However,  for  simplicity  as 
well  as  consistency  with  the  current  state-of-the-art  in  seakeeping  per¬ 
formance  assessment,  a  JONSWAP  expression  which  is  dependent  only  on  the 
two  parameters,  significant  wave  height  and  modal  wave  period,  is  also 
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desirable.  In  order  to  achieve  this,  the  following  steps  have  been 
carried  out.  Eliminating  the  dimensionless  frequency  f  from  equations 
(3)  and  (5), 
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The  fetch  dependence  of  the  dimensionless  total  energy  in  the  wave  spectrum 
is  given  as 
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(9) 


Furthermore,  another  fetch  empirical  relationship  is  well  represented  by 
the  following  relation 


1.6  (10 
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where  mQ  Is  the  spectral  moment  of  order  zero.  Significant  wave  height 
(Cw) 1/ 3  can  defined  as 
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and  substituting  into  equation  (9)  gives 
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Equation  (19)  can  be  further  rewritten  to  be  more  compatible  with  the 
usual  Bretschneider  spectral  density  formulation  as  well  as  most  ship 
response  calculation  procedures  by  converting  f  to  u>,  the  circular  frequency 
in  radians  •  sec  and  taking  y  ■  3-3, 


where 


or 


0  ■  0.07  for  <_  y~ 
o 

a  =■  0.09  for  >  f- 


(21) 

(22) 


B  PARAMETER 

The  modified  JONSWAP  spectrum  given  by  equation  (20)  contains  more 
energy  than  it  is  theoretically  supposed  to  for  the  same  values  of  sig¬ 
nificant  wave  height  and  modal  wave  period  when  these  two  parameters  are 
used  to  define  the  spectrum.  Figure  2  provides  an  illustration  of  the 
difference  in  significant  wave  height  at  the  same  fetch  for  winds  of  20 
and  30  knots.  The  solid  line  represents  the  theoretical  relationship 
between  significant  wave  height  and  fetch  for  those  wind  speeds.  The 
dashed  line  represents  the  values  which  are  actually  computed  from  the 
spectral  area  when  the  given  fetch  and  wind  speed  are  specified  in  the 
usual  JONSWAP  formulation  given  in  equation  (1).  The  difference  between 
the  solid  and  dashed  lines  represents  a  rather  noticeable  increase  in 
significant  wave  height  for  fetches  above  kO  nautical  miles.  The  diffi¬ 
culty  arises  due  to  the  fact  that  the  absolute  value  of  a  is  not  universal, 
as  was  first  suggested  by  Phillips,  and  also  as  described  by  Ewing  in 
Reference  k.  While  the  f  and  X  data  are  well  defined  in  the  linear  re¬ 
gression  fits,  see  equations  (5)  and  (10),  there  is  a  large  scatter  for 

£ 

This  inconsistency  is  also  present  when  the  usual  JONSWAP  formula¬ 
tion,  given  in  equation  (1),  Is  applied  and  is  particularly  noticeable  for 
fetches  above  about  kO  nautical  miles. 
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absolute  values  of  X  larger  than  10  ,  Consequently,  the  basic  form  of 
equation  (2),  and  hence  equation  (18),  is  not  universal  at  relatively  high 
waves  and  periods.  To  correct  this  anomaly,  a  new  parameter,  6,  is 
developed  to  replace  the  a  parameter  given  in  the  usual  JONSWAP  formulation. 
The  technique  developed  to  compute  $  Is  given  In  the  Appendix.  Figures  3 
and  4  are  the  result  and  provide  3  for  given  values  of  significant  wave 
height  and  modal  wave  period.  The  wave  parameter  ranges  are  deliberately 
broad  In  anticipation  of  any  extreme  occurrences  such  as  those  reported  in 
the  North  Sea.  The  3  parameter  assures  that  the  specified  (input)  sig¬ 
nificant  wave  height,  see  equation  (18),  corresponds  to  that  which  is 
calculated  from  the  Integration  of  the  resulting  spectral  ordinates.  The 
modified  JONSWAP  formulation  given  In  equation  (20)  can  now  be  rewritten 
as 


m^-sec  (23) 

where  3  is  a  constant  dependent  only  on  the  significant  wave  height, 

(c  )i/  ,  and  the  modal  wave  period,  T  . 
w  /3  o 

APPLICATIONS  OF  THE  MODIFIED  JONSWAP  SPECTRUM 
Figures  5  and  6  provide  sample  JONSWAP  (modified)  spectra  for  sig¬ 
nificant  wave  heights  of  2,  3,  5,  and  7  meters  and  a  range  of  modal  wave 
periods.  As  is  the  case  with  the  Bretschnei der  formulation,  equation  (23) 
can  be  applied  without  special  regard  to  fetch  or  wind  speed.  The  user, 
upon  selecting  the  values  for  significant  wave  height  and  modal  wave 
period  (e.g.,  from  historical  wave  statistics,  climatology,  etc.),  deter¬ 
mines  3  from  Figures  3  or  4  and  thence  the  spectrum  can  be  developed. 
Figures  3  and  4  contain  a  limited  but  realistic  range  of  modal  wave 
periods  for  given  significant  wave  heights  for  fetch  limited  geographies. 

Table  1  provides  comparisons  of  significant  wave  heights  for  different 
values  of  y.  By  substituting  a  set  of  significant  wave  height  values 
(e.g.,  from  Table  1)  In  the  abscissa  scale  of  Figures  3  and  4,  3  values 


can  be  obtained  for  various  values  of  y.  Since  a  range  of  fetch  or  wind 
speed  values  used  to  obtain  6  could  be  of  interest  in  some  seakeeping 
analyses  (e.g.,  for  specific  ocean  areas),  Figures  7  and  8,  developed 
from  computer  program  J0N1  as  detailed  In  the  Appendix,  provide  a  com¬ 
parison  with  corresponding  height  and  period  ranges. 

As  with  the  usual  JONSWAP  formulation,  the  modified  expression  given 
in  equation  (23)  is  for  long-crested  seas.  The  cosine  squared  spreading 
function  can  be  used  with  the  spectrum  but  Insufficient  data  on  the 
directionality  of  wave  systems  in  fetch- I Imi ted  waters  is  available  to 
verify  its  applicability. 

VALIDATION  OF  MODIFIED  FORMULATION 

The  general  shape  of  the  modified  JONSWAP  spectra  agree  well  with 
the  shape  of  the  measured  spectra  reported  for  Argus  Island  in  Reference 
5.  Several  such  comparisons  are  given  in  Figure  9.  Furthermore,  a  total 
of  four  spectra  in  the  North  Sea,  hindcast  by  the  U.S.  Navy's  Spectral 
Ocean  Wave  Model  (SOWM) ,  see  Reference  6,  were  randomly  selected  for 
comparison  with  the  modified  JONSWAP  formulation.  The  hindcasts  were  also 
compared  with  Bretschneider  spectra  defined  using  the  hindcast  significant 
wave  height  and  modal  period.  Figure  10  shows  these  comparisons.  In 
general,  the  modified  JONSWAP  spectra  provide  a  much  closer  approximation 
to  the  hindcasts  than  do  the  Bretschneider  spectra  though  secondary 
spectral  peaks  are,  of  course,  not  well  approximated  due  to  the  unimodal 
restriction  of  the  model.  The  modal  wave  period  has  been  used  as  the 
defining  parameter  of  the  spectra  in  this  comparison.  However,  since 
some  hindcast  spectra  contained  multiple  peaks,  average  wave  period 
(e.g.,  zero  crossing  period)  may  permit  a  better  shape  definition  of  the 
theoretical  wave  spectra. 


CONCLUDING  REMARKS 

An  expression  of  the  mean  JONSWAP  spectrum  dependent  on  the  signifi¬ 
cant  wave  height  and  modal  wave  period  is  derived.  However,  as  is  the 
case  with  the  usual  fetch-dependent  JONSWAP  spectrum,  it  contains  too  much 
energy  for  fetches  above  about  I»0  nautical  miles.  The  inconsistency 
arises  due  to  the  fact  that  the  absolute  value  of  a  is  not  universal. 
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while  the  f  and  X  data  are  well-defined  in  the  least  squares  fits 

(equations  (5)  and  (10)),  Further,  there  is  a  broad  scatter  for  absolute 

k 

values  of  X  larger  than  10  . 

B  is  developed  to  replace  the  a  parameter  given  in  the  usual  mean 
JONSWAP  formulation  and  to  correct  for  the  parameter’s  nonuniversality. 
Figures  3  and  4  permit  the  determination  of  B  for  given  values  of  signifi¬ 
cant  wave  height  and  modal  wave  period.  Hence,  wave  energy  is  conserved 
in  the  modified  JONSWAP  formulation. 

Recent  hindcast  spectra  from  the  North  Sea,  as  well  as  actual  wave 
measurements  from  fetch-limited  areas,  suggest  that  the  modified  JONSWAP 
spectrum  may  describe  wave  growth  conditions  more  realistically  than  the 
Bretschneider  spectrum  in  fetch-limited  or  shallow  water  conditions.  In 
general,  it  is  concluded  that  the  modified  JONSWAP  spectrum  provides  a 
reasonable  representation  of  wave  conditions  In  fetch-limited  or  shallow 
water  ocean  areas. 
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APPENDIX 

DESCRIPTION  OF  COMPUTER  PROGRAM  JON1 

It  has  been  shown  in  the  preceding  text  that  there  is  a  noticeable 
discrepancy  between  the  theoretical  and  actual  significant  wave  heights  for 
fetches  above  40  nautical  miles  when  the  JONSWAP  spectrum  is  applied. 
Therefore,  a  new  parameter,  0,  is  used  to  replace  the  a  parameter  in 
equatioh  (20)  >  eg, 


16.902  «w),/s 

1.375  T  2.75 
3  o 


1.375 


(24) 


where  (£w)i/  is  the  theoretical  significant  wave  height  in  meters. 

A  listing  of  program  J0N1  is  presented  in  Table  2.  The  purpose  of 
the  program  is  to  calculate  3,  modal  wave  period,  and  actual  significant 
wave  height  under  the  spectral  area.  Results  are  plotted  on  Figures  3  and 
4,  with  6  against  the  actual  significant  wave  height,  (c.  )i/  »  and  the 
modal  period,  T  .  For  example,  if  the  significant  wave  height  is  4.08  m, 
and  the  modal  period  is  8  sec,  then,  by  reading  across  the  intersection  at 
4.08  m  and  8  sec,  6  Is  seen  to  be  about  0.0135.  Substituting  B  value  into 
equation  (23)  and  the  modified  JONSWAP  spectral  can  be  generated  by  varying 
the  values  of  u>.  A  typical  output  from  program  J0N1  is  presented  in  Table 
3. 
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TABLE  2  -  COMPUTER  PROGRAM  JON) 
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TABLE  2  -  COMPUTER  PROGRAM  JONl  (Continued) 
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TABLE  3  -  OUTPUT  FROM  PROGRAM  JON1 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH  ,J 

NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 

THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM  I 

INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE  ^ 

THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 


